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Abstract 

The light distribution on the focal surface of spheric mirrors designed for 
telescopes in the lens-less Schmidt configuration is calculated analytically 
using geometrical optics. 

This analysis was motivated by considerations of the design the design of the 
AUGER fluorescence detector Its geometrical parameters are used in the 
examples. 

1. Introduction 

The lens-less Schimdt configuration employs a circular diaphragm at the 
center of curvature of the mirror. In this configuration the coma aberration 
is a second order effect and the dominant aberration is spherical. 
We are interested in calculating the light distribution on the focal surface of 
the spherical mirror. This information is relevant for an optimal design of 
the mirror and for use in detector simulation. 

In FigUJa schematic of the mirror is drawn. The adimensional normalized 
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Figure 1; Schematic of the mirror with the reference coordinate system. 
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coordinates zr = i/r = and pr = are introduced. 

The mirror is spherical with curvature radius Rm = 3400 mm and size 

3500 mm x 3500 mm . 

The hght entrance is hmited by a circular diaphragm located at the mirror 
center of curvature with a semiaperture Rd = 850 mm that defines an open- 
ing angle sin 6m = Rci/Rm- 

The light distribution on and close to the focal surface due to paraxial rays 
is calculated below including the effect of camera obscuration. 

2. Light distribution from paraxial rays without obscuration 

We consider a ray parallel to the zr axis at radial distance Hr = h/Ru < 
Rd/Ru- In Fig{T], the ray lies in the ur — zr plane. 

The ray intercepts the mirror at the point ((cos — |), sin B) where G is the 
angle with respect to the circle center (not the axis origin) . The reflected ray 
forms an angle 2B with the zr axis. Its equation in cylindrical coordinates 
PrZr is 



2.1. Calculation of the caustic 

The calculation of the light distribution on the focal surface requires the 
knowledge of the caustic of the spherical surface. The relevant formulae are 
in the Appendix. 

The relevance of the caustic stems from the fact that, for sin G < sin Gm, it is 



PR - sin G 



PR 
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the outer envelope of the converging rays. We are interested in studying the 
behavior of the rays close to the focal surface, where the previous condition 
on O is satisfied. Hence, the outer diameter of the light disk Pr{zr) is the 
maximum between the radial coordinates of the caustic and of the rays from 
the mirror rim, that is Eq{T]for O = Qm- 

2.2. Light envelope and disk of minimal confusion 

The size and the intensity of the light envelope are calculated by refor- 
mulating the equation of the caustic and solving the following system: 
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The solution is anticipated to be at |pr| << 1, therefore the approxima- 
tion in EqUHin the Appendix can be applied, resulting in the following third 
order polynomial equation: 
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Eqj3] can be solved with the Cardano method. 
The refiected ray from the largest angle 0m (mirror rim or diaphragm limited 
circle) intersects the caustic at two points. It is tangent at larger zr, by 
definition of the caustic, and has a simple intersection at smaller zr. At this 
point the radius of the light envelope is minimum. This is the disk of minimal 
confusion, given by 
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In Figl2] the envelope of light rays near the mirror focus is pictured to- 
gether with the caustic and a bundle of rays reflected from the mirror. The 
outer envelope is given by the rays converging from the outer rim for 

1 



2 V -1 + Vl - sin' 6m(1 + 2 sin' 6m)) < zn 



and by the caustic for 



ZRdm <zr< -(-1 + yl - sin'6Af(l + 2 sin' 6m)) 
and by the rays diverging from the outer rim for 

Zr < ZR4m 

where ZRdm is the position of the disk of minimal confusion obtained by Eql3l 
The radius of the outer envelope versus the radial distance is displayed in 
Figl3] (solid line) for Rm = 3400 mm. The envelope contains all the light 
reflected from the mirror. 

The position of the disk of minimal confusion is R = 1742.52 mm and its 
radius 7.32 mm. 
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In j2l-j3l, the AUGER design is analyzed using both user-written and pro- 
fessional ray tracing programs. The most straightforward comparison with 
these calculations is with where exactly the same geometrical parameters 
are used. The positions of the focal surface and the radius of the disk of 
minimal confusion agree within ^0.1 mm. 

2.3. Light distribution 

The radial light distribution at fixed zr is obtained by relating the 
known light distribution in to that in p^. The relation stems from 

df df dQ df dhji dQ df dA dhn dQ , . 

dpR dQ dpji dhji dQ dpR dA dh^ dQ dp^ 

The first factor is the inverse of the effective area of the mirror = 
{tiRm sin^ Qm)~^, the second is simply = 2TiR\jhji and the third is from 
= sin0, = COS0. The fourth is derived by differentiating EqU] 

dQ 1 

rln ~ 2 r ^ I In 1 ^^ tan 26 sine" " 

"'f^R cos2 2e'^^-R 2'^-^ cose'''' 2cos2e 
Substituting in Eq|5]we find 

df sin 26 1 



(6) 



dpR sin' 9m -^{zr + 1(1 - ^)) - 

rJi i« cos^ 2e ^ ^ 2\ cose^' 2cos'^e 

where 6 is expressed as a function of pr by inverting EqlH 

If more than one Q{pr) satisfies Eq{Tl the right side of Eql6] becomes the sum 

of all the solutions. 

The inversion requires the canonical change of variable t = tanG/2. In this 

and cos 6 = jj^. 



variable sin 6 = and cos 9 = tt^, so that EqU] becomes 



f + ^ " ' ' t' - - + 1 = (7) 

PR PR 
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Figure 2: Envelope of light rays near the focus without obscuration. Thick lines are 
external, dashed internal. 
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Figure 3: Spot radius versus radial distance from the focus: unobscured case (solid), 
obscured case (dashed). 
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This equation has 4 complex solutions and it is guaranteed that an even 
number of them are real although that might appear counterintuitive when 
Figl2] is analyzed graphically. 

In the region (marked with (3) in Figl2]) within the caustic and within the 
rays from the mirror rim after their crossing, three rays meet at a given p^?. 
One is converging at small G and has not yet reached the caustic, one is con- 
verging at larger G and has already reached the caustic and one is diverging 
at even larger and opposite sign G and has already reached the caustic on 
the opposite site. 

In the region (marked with (2) in Figl2]) within the caustic and outside the 
rays from the mirror rim, there are two rays, the diverging ray is missing. 
In the region (marked with Q) in FigE]) within the caustic and within the 
rays from the mirror rim before their crossing, there is only the converging 
ray that has not yet reached the caustic. In the region (marked with (T) in 
Figl2]) outside the caustic and within the rays from the mirror rim after their 
crossing, there is only the diverging ray that has already reached the caustic. 
In the region outside the caustic and outside the rays from the mirror rim, 
there are no rays. 

The reason for this apparent contradiction is that Eql7]for > always has 
one solution for t < — 1, that is G < — 7r/2, and vice versa for pji < 0. That 
corresponds to a ray refracted (not reflected) by an unphysical hemisphere 
specular image around the axis of the physical one. In other words, EqJT] 
represents full straight lines, while the reflected rays are only half straight 
lines. This solution is unphysical and must be neglected. 
Analysis of EqjT] leads to the conclusion that, for > 0, there is always an 
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additional real solution for — 1 < t < and a pair of solutions that can be 
either both real positive or complex conjugates. 

The reason for the presence of a region with two solutions is that the solu- 
tion at large and opposite sign 9 lays outside the physical region {—Qm, ©m)- 
The same reason requires the region with one solution of converging ray at 
small 0. 

In these two regions and in the region with three physical solutions, Eq|7] 
has four real solutions. In the region with one physical solution outside the 
caustic, EqI7]has two real and two complex solutions. 

EqJ7]can be solved exactly using the canonical Ferrari approach. Yet the solu- 
tion is very cumbersome and provides little insight into its physical meaning. 
The unphysical solution for the refracted ray can be removed since the ratio 
of the second to the fourth terms is where, at most, t = tanBM/2 ~ 
Gm/2 ~ 1/8. This ratio is at least 384 and dropping the fourth power term 
will affect the physical solutions only by a very small amount. With this 
approximation, EqlT] becomes 

The radial distribution ^ is plotted in FigHl The two-dimensional dis- 
tribution is plotted in FigjSl 

They are in good agreement with the light distributions shown in 

3. Light distribution from paraxial rays with obscuration 

In the AUGER design there is a 0.92 x 0.92 m^ photomultiplier camera 
with a spherical surface next to the mirror focal surface, that obscures a 
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Figure 4: Radial light distribution at the disk of least confusion for unobscured (solid) 
and obscured (dashed) cases 
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Figure 5: Two dimensional light distributions at — 10, 0, +10, +20 mm from the plane of 
least confusion 
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fraction of the incoming rays. 

The camera does not have rotational symmetry but its obscuration can be 
approximated by a disk with the same area, that is with radius Rc = 0.519 m 
positioned next to the focal surface zpt = 1/2. That implies that the rays hit 
the mirror only for angle sin 6 > sin 6c = Rc/Rm, that is for O > Oc = 8.78°. 

3.1. Light envelope and disk of minimal confusion 

The obscuration changes the light distribution along zr and might change 
the size and position of the disk of minimal confusion. The outer envelope 
of reflected light can be deduced by referring to FiglH] and noting that the 
light rays close to the zr axis are obscured. Hence the outer light envelope 
is given by the converging rays from the mirror rim up to the caustic, then 
by the caustic, then, possibly, by the converging rays from the rim of the 
obscuration disk, then by the diverging rays from the mirror rim. 
It depends on the obscuration area if the converging rays from the rim of the 
obscuration disk intersect the caustic at zjidQc) before or after zjidmi^Ai) 
from EqJH that is if 



Therefore if sin 6c < Y-2Mm(0A/)4/3, the outer envelope is unchanged, oth- 
erwise for ZRdm < zr < Zrc it is given by the converging rays from the rim of 
the obscuration disk. In this case the disk of minimal confusion is obtained 
by the intersection of the diverging ray from the outer rim with the converg- 
ing one on the opposite side from the rim of the obscuration disk. The result 
is 




sin^ Oc = zrc{Qc) 



where the caustic approximation of Eq{T4] is used. 
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This condition is verified in the AUGER design and the location of the 
disk of minimal confusion changes slightly. Its position is i? = 1742.80 mm 
and its radius is 7.17 mm. The radius of the envelope versus the radial dis- 
tance for the obscured case is displayed in FigJS] (dashed line). The difference 
between it and that of the unobscured case is minimal. 



3.2. Light distributions 

The results of the unobscured case are still valid, provided that only 
solutions with 6c < |6| < 6a/ are considered. 

The radial light distributions for the obscured case is presented in FigJH . 
4. Conclusions 

Detailed analytical calculations of the spherical aberration of a telescope 
spherical mirror designed for detecting fluorescence light emitted in the at- 
moshpere by extended air showers generated by very high energy cosmic ray 
interactions are presented. The position of the focal surface of the mirror 
and the light distribution on and close to it have been calculated. The cal- 
culations are useful both for the design of the detector and for the detector 
simulation. The AUGER Fluorescence Detector has been used as example 
but the results are applicable to any similar detector. 
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Figure 6: Envelope of light rays next to the focus with obscuration. Thick lines are 
external, dashed internal. 
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Appendix 

Third order equation 

We recall the solution of the third order polynomial equation. The original 
equation and the intermediate steps leading to the solution xi, X2 and x^ are 
listed (where u = v^, to ^ 1) 

= x^ + aix"^ + a2X + 03 
a? 2a? 0102 



P 



Q 



q q^ 



^2 V27 4 



Xi = P + g-^ X2 = ujP + uj^Q-^ x^ = uj^P + ujQ~'^ (10) 

3 2 

Two roots coincides if and only if |^ + y = 0. In this case, the third 
distinct root is Xi, that is the solution determining the disk of least confusion 
in our application. 

Some formulae about the caustic 

The calculation of the analytical form of the caustic from a spherical sur- 
face is a standard calculation, but it is surprisingly difficult to find it in text 
books. Following j4|, the caustic is the locus of intersection of neighbouring 
refiected rays when their distance becomes infinitesimal. Referring to the 
coordinate system in Eq{T] and using the dimensionless variables zr = 
and yji = it can be calculated by equating the derivative of Eq{T]with 
respect to O to and solving the system 
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The resulting curve can be expressed in parametric form as 



ynie) = sin^e 

zr{Q) = ^(-1 + \/l - sin2 6(1 + 2sin2 9)) (12) 

The caustic in three dimensions is obtained by rotating the curve in EqIT2] 
around the zr axis. 

It is useful to express exactly the caustic in alternative forms as zr^i/r), 
yR^ZR) and relating the parameter to ur and zr. From EqJT2l 



Q{yR) = aTcsm{^/y^) 



ZRiVR) = 2(-l + Vl-l/f (l + 2yf))- (13) 
A common and useful approximation for |sinG| ^ 1 (equivalent to 

I 1 2 

\yR\ <^ 1 and \zr\ <^ 1) is obtained from y 1 — yf^ ^ 1 — ^yf^, that implies 
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The exact expression for Ur^zji) or, equivalently, for Q{zji) is obtained 
from Eq{T2]by writing 

^^(9) = i(-l + cose(3-2cos2e)) 

which can be reformulated as a third order polynomial equation in cos G 

3 1 
cos^ Q - - cosQ + {zr + -) = 0. 

The discriminant of this equation is 

which is negative for — < zr < = zro. For these values, there are 
three real solutions: 



cos 



cos 



cos 



Q,{zr) = V2 



62 (^fi) = V2 



Qs{zr) = V2 



cos 



COS 



COS 



3n + arctan 




5n + arctan 


3 




TT + arctan ^ 


3 




3 

/ 



(15) 



The physical meaning of these solutions can be understood by considering 
them for Zr = 



cos0i(O) 
cos 62(0) 

cos 63(0) 
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The first is unphysical because the modulus is larger than one, the second 
corresponds to the tip of the caustic and the third is a large angle solution. 
Depending on the quadrant, the solutions are always monotonic versus zr 
either increasing (0i and 63) or decreasing (82) within the range — | < 
zr < zrq. The lower bound of the range is physical, the upper one comes 
from requiring three real solutions. For zr — zro, ©2 = ©3 = 7r/4 and for 
zr > Zro, ©2 and ©3 become complex conjugate. 

The solution relevant for the application to a mirror is therefore ©2 and we 
can write the exact parametrization of the caustic for < zr < zrq 
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The large angle solution is ©3 and reads 



Q{zr) 



VRizR) 



± arccos 



± 



— cos 2 



V 



-\/2cos 
/ 

V 



V 



. + arctan(^W^\\^ 



Acknowledgments 



We warmly acknowledge Prof. William Molzon from the Department of 
Physics and Astronomy of the University of California, Irvine, CA, USA, for 



19 



his careful proofreading of the manuscript. 
References 

[1] J. Abraham et al. Properties and performance of the prototype instru- 
ment for the Pierre Auger Observatory. Nucl. Instr. and Meth., A533:50- 
95, 2004. 

[2] G. Matthiae and P. Privitera. The Schmidt telescope with corrector 
plane. GAP-98-039. 

http://www.auger.org/admin/GAP_Notes/GAP1998/GAP_98_039.ps.gz. 

[3] M. Hrabovsky et al. The optical analysis of the proposed schmidt camera 
design. GAP-99-025. 

http://www.auger.org/admin/GAP_Notes/GAP1999/GAP_99_025.pdf. 
[4] B. Rossi. Ottica, chapter 2.7. Masson Italia Editori, Milano, 1984. 



20 



